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In order to explore the possible existence of the exotic state, we have constructed the 
tetraquark interpolating operators systematically. As a byproduct, we notice the 0" 1 tetraquark 
operators without derivatives do not exist. The special Lorentz structure of the 0~~ currents for- 
bids the four-quark type of corrections to the spectral density. Now the gluon condensates are the 
dominant power corrections. None of the seven interpolating currents supports a resonant signal. 
Therefore we conclude that the exotic state does not exist below 2 GeV, which is consistent 
with the current experimental observations. 
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I. INTRODUCTION 

Most of the experimentally observed hadrons can be interpreted as qq/qqq states and accommodated in the quark 
model Up to now there has accumulated some evidence of the exotic state with J PC = 1 + [1, 0, @ • Such a 

quantum number is not accessible for a pair of quark and anti-quark. It is sometimes labelled as an exotic hybrid meson 
with the particle contents qg s G^ v ^j v q. Recently we have investigated the 1 ^ state using the tetraquark currents 0). 
The extracted mass and characteristic decay pattern are quite similar to those expected for the exotic hybrid meson. 
Such a result is expected. Since the gluon field creates a pair of qq easily, the hybrid operator qg s G ^ v ^f q transforms 
into a tetraquark interpolating operator with the same exotic quantum number. In quantum field theory different 
operators with the same quantum number mix and tend to couple to the same physical state. 

Using the same tetraquark formalism developed in the study of the low- lying scalar mesons Q and the exotic f h 
mesons 0, we study the possible J PC = 0~~ states composed of light quarks. For a neutral quark model state qq, 
we know that J = ensures L = S hence C = {—) L+S = +1. In other words, states with J PC = 0~~, H are strictly 
forbidden. On the other hand, the gauge invariant scalar and pseudoscalar operators composed of a pair of the gluon 
field arc g*G« v G a > lv and e^^g^G^G^, both of which carry the even C-parity. 

We construct all the local tetraquark currents with J PC = 0~~. There are two kinds of constructions: (qqVqq) and 
(qq) (qq) . They can be related to each other by using the Fierz transformation. As usual, we use the first set [7fl . Their 
flavor structure can be 3/ (g) 3/, 6/ <g> 6/, and 3/®6/©6/<8>3/ ((qq)(qq)). With all these independent currents, we 
perform the QCD sum rule analysis. As a byproduct, we notice that there does not exist any tetraquark interpolating 
operator without derivative for the J PC = H case. 

This paper is organized as follows. In Sec. [TTJ, we construct the tetraquark currents with J PC = using the 
diquark (qq) and antidiquark (qq) fields. The tetraquark currents constructed with the quark-antiquark (qq) pairs are 
shown in AppcndixfA"! We present the spectral density in Sec. IHII and perform the numerical analysis in Sec. IIVI For 
comparison, we present the finite energy sum rule analysis in the Appcndix[B] The last section is a short summary. 

II. TETRAQUARK INTERPOLATING CURRENTS 

A. The J PC — 0~~ Tetraquark Interpolating Currents 

In this section, we construct the tetraquark interpolating currents with J PC = using diquark and antidiquark 
fields. Such a quantum number can not be accessed by a qq pair. The currents can be similarly constructed by using 
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the quark-antiquark pairs. However, as shown in Appendix^ these two constructions arc equivalent as we have 
shown several times in our previous studies @, 0] • 

The pseudoscalar tetraquark currents can be constructed using five independent diquark fields, which are con- 
structed by five independent 7-matrices 

Sabcd = {qia C( l2b)(q?,clbCql d ) , 

Vabcd = (qIa c i5q2b)(q3cCqJ d ) , 

Tabcd = (q^a Ca ^q2b){q?,ce^i<iCql d ) , (1) 
Aabcd = (qIa c i^q2b){q3ci' 1 i5CqJ d ) , 
Pabcd = {qI a C-f IJ .j5q2b)(q3cj fl CqJ d ). 

where qi_4 represents the up, down and strange quarks, and a — d are the color indices. 

To compose a color singlet pseudoscalar tetraquark current, the diquark and antidiquark should have the same 
color and spin symmetries. So the color structure of the tetraquark is either 6 © 6 or 3 <E> 3, which is denoted by labels 
6 and 3 respectively. Therefore, considering both the color and Lorentz structures, there arc altogether ten terms of 
products 



Wc list them as follows 



orentz 

®{3©6} 

Color ■ (2) 

q^a C q2b(q3albCqJ b + q3blbCqJ a ) , 
6 F <g> 6 F (S) { Vq = qla C J5q2b(q~3aCqJ b + q~3bCqJ a ) , 

qL C,J ^q2b{q3aCr^75CqJ b - q3b^ v ^f5CqJ a ) , 

qTa C q2b{q3alsCql b - q 3 blbCqJ a ) , 
3 F ® 3 F (A) { V 3 = q{ a Cj5q2b(q3aCqJ b - q~3bCqJ a ) , (3) 

qI a C(7p. V q2b{q3a0^1bCql b + ^b^lbCqJa) , 
qL C l^q2b(q3a^l5CqJ b + g 3 67 M 75CgJ a ) , 

qTaCinJ5q2b(q3aj t *CqJ b - q 3 bY l CqJ a ) . 

6 F ®3 F (M) { 7lV 



3 F ® 6 F (M) p 6 _ _ (s _ ^ nsT s ^ n # 



qiaCl f ,l5q2b(q3al> J 'Cq 4b + g 3 b7 M Cg 4a ) , 
qia C l^Q.2b{q.3al^lhCqJ b - qsbl^lhCqJa) 



In the above expressions, the flavor structure is fixed at the same time due to the Pauli principle. The currents Sq, 
Vq, T 3 belong to the symmetric flavor representation 6p <8> 6f(S) where both diquark and antidiquark fields have the 
symmetric flavor structure. The currents S3, V 3 , Tq belong to the antisymmetric flavor representation 3f <8> 3f(A), 
where both diquark and antidiquark fields have the antisymmetric flavor structure. Aq, P 3 for 3f(S>6f(M) and A3, P@ 
for 6p ® 3p(M), where M represents the mixed flavor symmetry. The isovector states with charges can be observed 
in the experiments more easily, therefore in thispaper we concentrate on the isovector currents which was shown in 
the S/7(3) tetraquark weight diagram in Fig 1 Q. We have: 

qqqq(S), qsqs(S) € 6 F ® 6 F (S) , 

qsqs(A) e 3 F ®3 F (A), (4) 
qqqq(M),qsqs(M) € (3 F ® 6 F ) © (6 F ® 3 F ) (M ) . 

We do not differentiate up and down quarks and denote them by q. We are only interested in those neutral components. 
The other states do not carry definite C-parity. It turns out that the neutral isovector and isoscalar states have the same 
QCD sum rules. Our following discussions are valid for both of them. Taking the charge- conjugation transformation, 
we get 

cs 6 cr 1 = Vq , CAgCr 1 = p 6 , ca 3 c _1 = p 3 , cs 3 <cr 1 = v 3 , CT 6 cr 1 = r 6 , CTgtcr 1 = r 3 . (5) 

Tq and T 3 have even charge-conjugation parity. We conclude that the currents with J PC = arc: 

l!a C q2b(q3al5CqJ b + q~3blhCql a ) - qia C J5q2b(q~3aCqJ b + q~3bCqJ a ) , 





= s 6 


-V e 


r ] 


= A 6 


-Pe 




= A 3 


-P3 




= s 3 


-v 3 



qIaCl f ,q2b(q3al' Jl l5CqJ b - 9367^75 CqJ a ) - qlaClnlsq2b(q3aY L Cql b - q3bl lL Cq T Sta 
lIa C q2b(q3al5CqJ b - q~3blhCql a ) - qia C lzq2b{q3aCql b - q~3bCqJ a ) ■ 



(6) 



Adding different quauk contents as shown in Eq. ([¥]), there arc altogether seven independent currents as shown : 
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1. For 6 F ®6 F (S): 

??i = S 6 (qqqq) - V 6 (qqqq) = u^Cd b (u a j 5 CdJ + u b "f 5 Cd%) - u^Cj 5 d b (u a CdJ + u b Cd^) , 

V2 = S 6 (qsqs) -V 6 (qsqs) = u^Cs b {u a ^ 5 Csl + %y 5 Cs^) - u^Cj 5 s b (u a Csl + u b Cs^) , (7) 

where 771 belongs to the 27p representation and contains up and down quarks only while 772 belongs to the 8f 
representation and contains one ss quark pair. 

2. For (3 F (g> 6 F ) © (6f <S> 3 F ) (Af): 

% = A 6 (qqqq) - P 6 (qqqq) = u^C^^d b {u a ^^CdJ + u b ^^ 5 CdJ) - u^C^^ b d b {u a ^CdJ + u b ^CdJ) , 

m = A 6 (qsqs) - Pe(qsqs) = vF a Cj^s^Ua^^Csl + u^^Cs^) - u^C^^s^u^Csl + u^Cs^) , 

775 = A 3 (qqqq) - P 3 (qqqq) = ulC^d b (u a ^Cd[ - U b -fC%) - u T a C ltll5 d b {u a ^CdJ - Ub^Cd^) , (8) 

V6 = A 3 (qsqs) - P 3 (qsqs) = u^C-f^s b {u a YCsJ - u b ^Csl) - u^Cj^ 5 s b (u a 'y' J 'Cs b - u^Cs^) , 

where 773 and 775 belong to the 10f representation and contain only u, d quarks while 774 and r]e belong to the 
8f representation and contain one ss quark pair. 

3. For 3 F ® 3 F {A): 

V7 = S 3 (qsqs) - V 3 (qsqs) = u^Cs b (^75 Cs b - u b ^Cs^) - u^C^ 5 s b {u a Csl - u&Cs^ ) • (9) 
where 777 belongs to the 8f and contains one ss quark pair. 
It is understood that there exists the other part ±[u <->■ d] in the expressions of 772,4,6,7- 

B. The J PC = +_ Tetraquark Currents 

Now we move on to the J PC = 0^ case. There are also ten independent scalar tetraquark currents without 
derivative: 

S 6 = Qia C( l2b{q3aCqJ b + q 3 bCqJ a ) , 

V 6 = qial^C^biqzaC^qJi, + qsbC^ql a ) , 

^6 = qL<J^Cq 2 b(q 3 aCa^qJ b + q 3b Ca^qJ a ) , 

A'e = qIal^l5Cq2b(q 3a C^J 5 qJ b + q~3bC^J 5 qJ a ) , 

P 6 = Qial5Cq2b(q3aC'y 5 qJ b + q 3 bC^ 5 qJ a ) , 

S 3 = qIa C( l2b(q3aCqJ b - q 3b CqJ a ) , (10) 
V3 = qlalnCq 2b {q 3 aC^ql b - q 3 bC^qJ a ) , 
T-3 = qL^Cq 2b (q 3a Ca^qJ b - q 3b Ca^qJ a ) , 
A3 = qIalnl5Cq 2b {q3aC^J 5 qJ b - q~3bC^J 5 qJ a ) , 
P 3 = qiaJ5Cq 2 b(q3aCj 5 qJ b - q~3bC-f 5 qJ a ) . 

The flavor structure is again fixed due to the Pauli principle. To have a charge-conjugation parity, we fix the quark 
contents to be: q\ = q3 and q 2 = (74 (or q\ = 54 and q 2 = (73). After performing the charge-conjugation transformation, 
we find that they all have an even charge- conjugation parity, for example: 

CS'qC- 1 = +S' 6 . (11) 

Therefore, the J PC = 0^ tetraquark interpolating currents without derivatives do not exist. 
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III. THE SPECTRAL DENSITY 

We consider the two-point correlation function in the framework of QCD sum rule [1, Q: 

n(q 2 ) = f d 4 xe lqx (Q\T V (x)rf(0)\0), (12) 



where r\ is an interpolating current. We can calculate II(g 2 ) at the quark gluon level using the propagator: 

iSf ee <0|T[g»g- & (0)]|0> 

i5 ab . i A", _ 1 , S ab 



2tt 2 x 4 32. 2 2 * ^x 2 " ' 12 

b abr r 2 m q S ab i5 ab m q {qq) . zg o6 m^gg) . 

(g s gCTGg) - — ^ + x + ————x, (13) 



192 ^ ^ Hl 4tt 2 x 2 48 8tt 2 x 

where x = juX**. With the dispersion relation n(g 2 ) is related to the observable at the hadron level 

n(p 2 ) = r p( 2 s) . ds, (14) 

Jo s~p 2 -is 

where 

p(s) ee ^ 5(s - M 2 )(0\r]\n)(n\rf\0) = f%5{s - M%) + continuum . (15) 

n 

Here, the usual pole plus continuum paramctrization of the hadronic spectral density is adopted. Up to dimension 
12, the spectral density Pz(s) at the quark and gluon level reads: 

/ \ ™\ 3 ( (g 2 GG) m q (gg) 2 (g 2 s GG)m 2 q (tffGGG) s 

P1[S) ~ 15360^ "192^ S " ( 3072^ ~ ~2A^ )& +[ 256.6 + 768.6 (^(^2) ~ 5)J* 

3m 2 (qq} 2 {g 2 GG)m q (qq) 16 3 1 2 _ 

~^ — 2^2 — + 192^4 ) + y-^ m im) ~ ■^ m q m)(9sqoGq))5{s) , (16) 

/ n = g4 m s 3 , / ™ 4 , m s(ss) (g 2 GG) 2 

P2U 15360. 6 384tt 6 l 64. 6 48. 4 3072. 6 j * 

, r (^/GGG)_ (31n(A) _ 5) _ ( mf(S*> m 2 (g 2 GG)^„ 



768tt 6 v V y ' 8. 4 512. 6 
mg(ss) 2 m 2 s {qq) 2 m s {ss)(g 2 GG) m 2 {uu){g B qaGq) 8 2 
+ ( ^ 3^ 384^ } - ( 6^ 9»».<">(W) > ( 17 ) 



M * 4 m 2 q 3 5{,g 2 GG) m,(gg) 2 AglfGGG) s b(g 2 s GG)m 2 

P3{S) = 3840^~48^ S + ( T536^ + ^^ )S +[ 192.6 ( 3hl (^)-5)- 12Bn * I s 

, 6mg<gg) 2 5(g 2 GG)m q (qq) , _ ,64 3 4 2 _ 

-( ^2 T^4 ) + (y m 9<W - ^rn q (qq)(g s qaGq))5{s) , (18) 

/ \ s 4 3 m 4 TO s (sa) 5(.g 2 GG) 2 

P4l ' J 3840tt 6 96. 6 l 16. 6 12. 4 1536. 6 j " 

MPJGGG) s m*(s S ) 5m 2 s (g 2 GG) 

+ [ 192^ (31n( ^ } - 5) - + 256^ )]S 

, ( m s(ss) 2 _ 4m 2 (gg) 2 5m s (5 S )(g 2 GG) _ 2m 2 (^)(g,ggGg) _ 32 / 5 .\/ ff „\2 U /^ / 1q n 
+l 3. 2 3.2 + 192. 4 j 1 3.2 9 »».WW«MaW, U»; 



2 3 ,(.g 2 GG) 777, (gg). 2 r (9sfGGG) . « . (5 2 GG)m 2 



7680. 6 96. 6 v 1536. 6 12. 4 y 384. 6 v ' 128. 6 

,2 



^(gg) 2 (g 2 GG)m q (gq) ^ , ,32 3 2 

( ^5 gg^x ) + - ^2m q (qq)(g s qaGq))5(s) , (20) 
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Pb(s) = 



7680tt 6 192tt 6 
, AglfGGG) 



* 3 + ( 



32tt 6 



+(' 



384tt 6 
6tt 2 



(31n(^)- 
M 

2m 2 s (qq) 2 
3tt 2 



5)-( 



m s {ss) 

24tt 4 

3 (ss) 



III 



(g 2 s GG) 



1536tt 6 



4tt 4 
^)(g s 2 GG) 

1927T 4 



)-( 



256tt 6 /J 

m 2 s (uu)(g s qaGq) 



3tt 2 



16 



«i s (ss)(<?g) 2 )<5(s) : 



(21) 



s 4 m 2 3 m 4 m s (ss) (g 2 s GG) 2 

Pn J 30720tt 6 768tt 6 l 1287r 6 96tt 4 3072tt 6 j 

+ [ 1536.6 (31n( ^ } - 5) ( ^6^ + 512^ )]S 

, f m 2 s (ss) 2 m 2 s {qq) 2 m s {ss)(g 2 GG) m 2 s {uu){g s qaGq) 4 2 
+ ( ^ 6^ + —§84^ } ( 12^ 9»».<">(M) ■ ( 22 ) 

It is interesting to note several important features of the above spectral densities: 

• First the special Lorentz structure of the J = interpolating currents forbids the appearance of the four- 
quark type of condensates (qq) 2 , (qq) (g s qaGq) and (g s qcrGq) 2 . Usually these terms play an important role in 
the multiquark sum rules. The Feynman diagrams for the dimension 10 condensate (g s qaGq) 2 are shown in Fig. 
ffl 




FIG. 1: Feynman diagrams for the quark gluon mixed condensate. 



• The dominant non-perturbative correction arises from the gluon condensate, which is destructive for pi-2(s) 
and constructive for P3_7(s). Moreover there are corrections from the tri-gluon condensate (g 3 f abc G a G b G c ) as 
shown in FigO In the above expressions we use the short-hand notation (g^fGGG) to denote the tri-gluon 
condensate. There are three types of Feynman diagrams. The first class of Feynman diagrams vanishes because 
of the product of the color matrices. The second class is proportional to m q and could be omitted in the chiral 
limit. Only the third class leads to non-vanishing tri-gluon correction. In fact the gluon condensates become 
the only power corrections in the chiral limit. 




FIG. 2: Feynman diagrams for the tri-gluon condensate. 



• The second term in each Pi{s) is destructive, which renders the spectral density negative when s is small. 
This —m 2 s 3 piece is an artefact of the expansion of the quark propagator ^_ 8 m in terms of the quark mass 
TTLq perturbatively. Without making such an expansion, the perturbativc contribution to the spectral density 
is always positive-definite. Such a destructive term will sometimes produce an artificial plateau and stability 
window in the sum rule analysis, which must be removed. 
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• Although the tree-level four-quark condensate vanishes, one may wonder whether the four-quark condensate 
g 2 s (qq) 2 plays a role since the latter is very important in the qq meson sum rules [H, Two types of Feynman 
diagrams could produce such a correction. The first class of Feynman diagrams is very similar to that in the qq 
meson case where a gluon propagator is attached between two-quark condensates, as FigJ3] shown. It's easy to 
check that they vanish due to the special Lorentz structure of the correlation function. One of the second class 




FIG. 3: One set of Feynman diagrams for the four-quark condensate. 

of diagrams is shown in Fig[¥] In this case, we use the mesonic type interpolating currents in the appcndbjAl to 
simplify the derivation. After making Wick-contraction to the correlation function, 

we get 

Tr[-r[SQ(x - y)r' 2 S Q (y - x)]Tr[-S Q (x - 22)7" ' Sq{z 2 ~ y)T 2 S Q (y - z 1 )^Sq(z 1 - x)Ti x g^ v x S G (z2 - Zi)]. 

where Sq is the quark propagator and Sq is the gluon propagator. {Ti^} could be cither {/, 75} or {j a , 
757a}- Sq(i/ — zi) oc (qq). In fact, there would be three 7-matrices or three 7-matrices plus 75 left in the latter 
trace. Therefore this piece also vanishes. 




FIG. 4: Feynman diagrams for the four-quark condensate. 



IV. NUMERICAL ANALYSIS 



In the chiral limit (m s 



P3-4(s) 

P5-e{s) 
Pr(s) 



0) the spectral density reads 

(g 2 sGG) g2 



15360tt 6 3072tt 6 



3840tt 6 
7680tt 6 

4 



(glfGGG) 
768tt 6 

H9 2 S GG) a , (glfGGG) 
1536tt 6 S 

(g 2 GG) 2 

1536tt 6 

(g 2 s GG) 2 



30720tt 6 3072tt 6 



192tt 6 
(glfGGG) 
384tt 6 
(glfGGG) 
1536tt 6 



(31n(— )-5)s, 
M 

(31n(4)-5)s, 
M 

(31n(4)-5)s, 
M 

(31n(^)-5) S 



(23) 



where jl = 1 GeV. Requiring the pole contribution is larger than 40%, one gets the upper bound M^^x of the Borel 
parameter Af|. The convergence of the operator expansion product leads to the lower bound M 2 • of the Borel 
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\ 


\M 2 ■ M?n a v~\(SVZ) 
I rnin ' ivJ rnaxj v° v ^) 


\M 2 ■ Mmavl Qoffel 
L rnin ' JlJ maxj ^lun^y 


Pl-2 


0.77 ~ 1.50 


0.90 ~ 1.68 


P3-4 


1.22 ~ 1.90 


1.40 ~ 1.65 


P5-6 


1.05 ~ 1.77 


1.55 ~ 1.74 


pi 


1.10 ~ 1.85 


1.50 ~ 1.75 



TABLE I: The working region of M% with Ioffe's and SVZ's ghion condensates and so = 7 GeV 2 . 



parameter. In the present case, we require that the two gluon condensate correction be less than one third of the 
perturbative term and the tri-gluon condensate correction less than one third of the gluon condensate correction. The 
working region of Mg in the sum rule analysis is [-^ mmJ M^ ax ], which is dependent on the threshold sq. 

In order to study the sensitivity of the sum rule to the condensate values, we adopt two sets of the gluon condensate 
values in our numerical analysis. One set is from Ioffe's recent review [13 : (g 2 s GG) = (0.20±0.16) GeV 4 , (g 3 JGGG) = 
0.12 GeV 6 . We also use the original SVZ values [|: (g 2 s GG) = (0.48 ± 0.14) GeV 4 , {g 3 JGGG) = 0.045 GeV 6 . The 
working regions of the sum rules with the above two sets of gluon condensates and so = 7 GeV 2 are listed in Table 
HI The working region of the sum rule is very narrow even with Sq = 7 GeV 2 . The variation of Mx with Mg and sn 
is shown in Figs. O5]for the interpolating currents 771—2 , 773-4, 775—6, respectively using Ioffe's gluon condensate 
values. The variation of Mx with M% and sq and SVZ's gluon condensate values is presented in Figsl9lfl~2l 

For a genuine hadron state, one expects that the extracted mass from the sum rule analysis is stable with the 
reasonable variation of the Borel parameter and the continuum threshold. In other words, there should exists dual 
stability in M| and sq in the working region of M|. From all these figures we notice none of the mass curves satisfy 
the stability requirement. These interpolating currents do not support a low- lying resonant signal. 
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FIG. 5: The variation of Mx with M% (Left) and so (Right) for the current 771-2 using Ioffe's gluon condensate values. 
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FIG. 9: The variation of Mx with Mg (Left) and so (Right) for the current 771-2 using SVZ's gluon condensate values. 
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FIG. 10: The variation of Mx with M% (Left) and so (Right) for the current 773-4 using SVZ's gluon condensate values. 
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FIG. 11: The variation of Mx with M B (Left) and so (Right) for the current 775-6 using SVZ's gluon condensate values. 
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FIG. 12: The variation of Mx with M.% (Left) and so (Right) for the current 777 using SVZ's gluon condensate values. 



V. CONCLUSION 



The exotic state with J = can not be composed of a pair of gluons nor qq. In order to explore the possible 
existence of these interesting states, we first construct the tetraquark type interpolating operators systematically. 
As a byproduct, we notice that the J PC = 0" 1 tetraquark operators without derivatives do not exist. Then we 
make the operator product expansion and extract the spectral density. The gluon condensate becomes the dominant 
power correction. Usually the four-quark type of condensates {qq) 2 , (qq) (g s qaGq) and (g s qaGq) 2 are the dominant 
nonperturbative corrections in the multiquark sum rules. However these terms vanish because of the special Lorentz 
structure imposed by the exotic quantum numbers. 

Within the framework of the SVZ sum rule, we note that the absence of the (qq) 2 , (qq) (g s qaGq) and (g s qaGq) 2 
terms destabilize the sum rule. There does not exist stability in either M% or sq in the working region of Mg. 
Therefore wc conclude that none of these independent interpolating currents support a resonant signal below 2 GeV, 
which is consistent with the current experimental measurement [l[. 
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APPENDIX A: INTERPOLATING CURRENTS IN (qq)(qq) BASIS 

For 6 F ® 6> (S): 

1m )(1) = {qial t j.qia){q2bl' 1 lbq2b) + (qiaj^5qia)(q2bj' 1 q2b) + {qial^aji&bl^lbqib) + (qiaj^5q2a)(q2bj ll qib), 

7 ?m )(8) = ^abK.d{(qiai^qib)(q2ci' J 'i5q2d) + {q\ai,iihqib){q2di' 1 q2d) + {qiai^bji&d^isqid) + {qiai^ibq2b){q2cY L qid)}, 

For (3f ® 6 F ) © (6 F ® 3 F ) (M): 

Vim = (qiaqia)(q~2b75q2b) ~ (9la753la)(?2&g , 2&) , 

^(M)(8) = X ab \ cd {(q la q lb )(q 2cl5q2d ) - (q la "f 5 q lb )(q 2c q2d)}, 

V2m = {qialnqia){q~2bl ll l5q2b) ~ (<Zla7 V^Ola) (<Z2b7 P ''12b) , 

l2m <8) = ^abKd{(qiaJnqib)(q~2cj' 1 j5q2d) - (q~la7^5qib){q~2c^q2d)}, 

For 3 F ® 3 F (A): 

Vm ){1) = (qiaj^qia)(q2bj ll j5q2b) + (qiaj^qia) {q2bi' 1 q2b) - (qiaj^q2a)(q2bj ll j5qib) - (qiaj^j5q2a)(q2bj ll qib), 

Vm = KbKd{{qiai tJ ,qib)(q2ci tl ibq2d) + (qiai^5qib)(q2dY l q2d) - {qiai^bji&ci^isqid) - {qiai^i^bji&ci^qid)}, 

where the indices (1), (8) represent the color singlet and octet. Now we get eight mesonic currents. Then we introduce 
the formula of the interchange of the color indices: 

(9i g , 2&93o?4&) = 7^{qiaq2bq~3bq4a) + ^ Kb^cd{qi a q2cq3dqib) , 
16 1 

^ab^cd(qia,q2cq3bq4d) = ~g (<Z la 92b <?3 b ?4a) — -^ab^cd(qiaQ2cQ3dq4b), (Al) 

Next, we perform the Fierz rearrangement in the Lorrentz indices with the formula 

(ab)(ba) = -(aa)(bb) + ^{a l5 a)(b l5 b) + ^{a^a){h"b) - ^{a l5 ^a)(b l5 ^b) + ^{a<j^a){b<j^b) 7 (A2) 
For example, we have 

(qI a Cq2b)(q3al5CqJ b ) = ~^{qi a C-i 5 Cql b ){q 3a q 2 b) - ^{ql a Cl^l5Cql b ){q3al^q2b) 

\(qIaC / ^^l5CqJ b )(q3aCr t " U q2b) + j(qLCl t iJ5'~f5CqJ b ){q3aj' 1 l5q2b) 



4 

" 1 (q[a C l5l5 C qj b ) (<73a75<?2fc) 
1 

t 



= -^(94b75<7la)(93a926) ~ ^(g4b7p759la)(93a7' i 92 ) (A3) 
+ o(94b^75<7la)(<Z3aCr' 1I/ <72fc) ~ J (^bl^la) (?3a7 A '7592b) 



-^(94b9la)(<73a7592b), 
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There are only four independent currents among those eight mesonic currents. Any four currents are independent 
and can be expressed by the other four. 



(S)(8) _ 



V 



'llm 

JM)(8) 
'hm 



,(«)(!) 



- 3 Vlm 

. 4 (A/)(D 
/lm 



(M)(l) 



(M)(l) 



(A)(1) 



We establish the relations between the diquark currents and the mesonic currents using the Fierz transformation. For 
instance, we can verify the relations 



n (S)(l) 
'llm 

'hm 
„(A)(1) 



\W l 4- -n M 
, 'lid + 4 '2d' 

_Iti m 4- -n M 
2 "id + 



APPENDIX B: FINITE ENERGY SUM RULE 



Sometimes the finite energy sum rule is also employed in the numerical analysis. One first defines the nth moment 
using the spectral density 



(Bl) 



W(n,s ) = / p{s)s n ds. 
Jo 

With the quark-hadron duality we have 

W(n, S )\Hadron = W(n, S )\oPE- 

The mass of the ground state can be obtained as 



M%(n,s ) 



W(n+l,3o) 
W(n,s ) 



(B2) 



(B3) 



We have plotted the variation of Mx with So for all the seven interpolating currents in Fig. [T3l The left and right 
diagrams correspond to Ioffe's and SVZ's gluon condensate values respectively. It seems that there exists a minimum 
of Mx for each current. However, a reasonable sum rule requires that the operator product expansion should converge 
well. In other words, we require that the two-gluon power correction be less than one third of the perturbative term 
and the tri-gluon power correction less than one third of two-gluon power correction in W(0, So), which leads to the 
working window of this finite energy sum rule as: 



\ 


so(SVZ) 


so(Ioffe) 


Pl-2 


4.0 


7.0 


P3-4 


4.2 


5.7 


P5-6 


4.0 


7.0 


P7 


4.9 


6.0 



Clearly for each current the minimum of the mass curve lies outside of the working region in both of the figures and 
is not a real resonant signal. Starting from 4.0 GeV 2 , each mass curve grows monotonically with sq. Thus, there does 
not exist a resonant signal for every interpolating current. 
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FIG. 13: The variation of Mx with so and n — from the finite energy sum rule. The left and right diagrams correspond to 
Ioffe's and SVZ's gluon condensate values respectively. 



